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Effects of Residual Stress on Static and Dynamic Characteristics of an 
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Electrostatically actuated nanotubes and nanowires have many promising applications as nano- 
switches, ultra sensitive sensors and signal processing elements. These devices can be modelled as 
slender beams with circular cross-section. In this paper, effects of residual stress on static and dy- 
namic characteristics of a cylindrical nanobeam are presented. Galerkin based multi-modal reduced 
order model technique has been used to solve the governing differential equation. The equation 
has also been solved numerically by collocation method for verification of the results. The effects 
of higher modes in reduced order modelling have been investigated. The analysis shows that the 
residual stress significantly influence the static and dynamic characteristics of the nanobeam. In 
particular, the results show that resonating frequency tunability increases significantly under the 
presence of compressive residual stress. In addition, an inherent resonating frequency stability 
mechanism under temperature variation in nanobeam resonators has also been explored. 
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I. INTRODUCTION 

Nanoelectromechanical systems (NEMS) is a vibrant 
field of research. Many NEMS devices have been de- 
veloped as ultrasensitive sensors [HjH, nano-switches @ 
and RF communication elements [4|, based on electro- 
statically actuated nanotubes and nanowires. In electro- 
static actuation, a nanotube or nanowire is in capacitive 
arrangement with a plate electrode and bias voltage be- 
tween them creates actuation force. DC voltage is re- 
sponsible for static deflection while AC voltage creates 
oscillatory motion. Efforts have been going on to un- 
derstand the various static and dynamic characteristics 
of nanotubes and nanowires based devices 0-0] ■ NEMS 
is a natural advancement of microelectromechanical sys- 
tems (MEMS). In NEMS, nanotubes and nanowires with 
circular cross-sections form a new class of beam struc- 
tures while in MEMS rectangular cross-section beams are 
prevalent. Many efforts have been devoted to understand 
the characteristics of microbeams 

Galerkin based reduced order modelling is an impor- 
tant technique which is widely used in MEMS analysis 
[1, El- In this method, the governing integro-partial dif- 
ferential equation (IPDE) of motion is reduced to a sys- 
tem of ordinary differential equations (ODEs). This is 
done by assuming the solution in terms of pre-assumed 
spatial functions (mainly linear mode shapes of a beam) 
and unknown time functions, and then substituting the 
assumed solution into the IPDE to obtain the required 
system of ODEs. In recent works, Ouakad and Younis 
and Rasekh et al. [lj| have carried out static and 
dynamic analysis of nanotube oscillators using reduced 
order model (ROM) technique. Rasekh et al. have used 
single-mode ROM for analysis while Ouakad and You- 
nis have used multi-modal ROM, and observed that only 



single basis function (first linear mode) is sufficient to 
obtain a converged solution. 

Residual stress, often introduced during fabrication, 
alters the designed characteristics of actuators and oscil- 
lators (l3j . Fixed- fixed type beams are dominated by ax- 
ial thermo-mechanical residual stress, originating in the 
difference in thermal expansion coefficients of the beam 
and supporting structure. Residual stress has been con- 
sidered by some researchers on rec tang ular cross-section 
microbeams. Abdel-Rahman et al. [14J have solved static 
and eigenvalue problem of microbeams using shooting 
technique. They included the effect of residual stress in 
their study. In eigenvalue analysis, free vibration charac- 
teristics of a microbeam at deflected position due to ap- 
plied DC volta ge h as been studied. Kuang and Chen [15[ 
and Jia et al. [16j | have done eigenvalue analysis of mi- 
crobeams, including the effect of residual stress, using dif- 
ferential quadrature method. Soma et al. have stud- 
ied experimentally and theoretically the effects of resid- 
ual stress in pull-in behaviour of micro-switches. Pull-in 
is a well-known instability point beyond which elastic 
restoring force cannot balance applied electrostatic force 
and beam snaps down to the gate electrode @. 

The focus of this study is a circular cross-section 
nanobeam with fixed-fixed end condition. The effects of 
both tensile and compressive type axial residual stresses 
have been studied using multi-modal ROM technique. 
The study comprises static deflection under applied DC 
voltage, pull-in parameters and variation of resonating 
frequency with applied DC voltage. The convergence of 
the ROM solution has also been studied. The results 
have been verified by numerically solving the boundary 
value problem using collocation technique. In addition, a 
mechanism is suggested for resonating frequency stability 
under temperature variation. 
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II. MATHEMATICAL MODEL 

A schematic diagram of an electrostatically actuated 
cylindrical nanobeam is depicted in Fig. [TJ The beam has 
clamped-clamped end condition with length L and radius 
R. It is in capacitive arrangement with a electrode plate 
where g is the distance of separation between beam and 
plate. The hat represents physical quantities with dimen- 
sion to differentiate from their non-dimensional forms. 
The beam is actuated by a DC load Vdc superimposed 
an AC harmonic load of amplitude Vac and driving fre- 
quency uif. The beam is vibrated in viscous medium 
which has the damping coefficient c. The residual stress 
is assumed uniform throughout the beam and is axial in 
nature. Here axial load N represents the residual stress of 
the beam. u(x,t) is vertical deflection of the nanobeam 
along spatial coordinate x at time t. The equation of 



where T = y pAL 4 /EI is a time constant. By substi- 
tuting the non-dimensional variables in (fTJ), the modified 
equation can be written as 



u + u + tit 
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The overhead dot ■ and prime represent partial deriva- 
tives with respect to time t and spatial coordinate x re- 
spectively. In the modified equation ([5]), AC and DC 
voltage terms are kept in dimensional form. Static anal- 
ysis of a nanobeam under the action of applied DC load 
Vdc is presented in next section. 



FIG. 1. Schematic diagram of an electrostatically actuated 
nanobeam 



III. STATIC ANALYSIS 
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and the boundary conditions are 



(1) 



du 

dx 



du 
dx 



\x= 



\x= 1 



= 0- 



A and I represent the cross-section area and the area mo- 
ment of inertia of the nanobeam. E and p are Young's 
modulus and density of the nanobeam and £o is the per- 
mittivity of vacuum. All terms in the left hand side of 
([T]) are same as in a equation of free vibrating linear 
beam. The first term of the right hand side of (fTJ) in- 
cludes nonlinear stretching term due to deflection and 
residual stress whereas the second nonlinear term is ac- 
tuation force due to applied DC and AC voltage. 

For convenience, ([T]) is represented in terms of non- 
dimensional variables 

u x ~ t 

u = — , x = — and t = — , 
9 L T 



Static response of a nanobeam whose properties is 
given in the Table U is studied. This nanobeam model 
has also been analysed by Ouakad and Younis Q and Ke 
et al. [l8j |. The equation for the static analysis can be 



TABLE I. Properties of the nanobeam which is used for anal- 
ysis 



Length L 
(nm) 


Radius R 
(nm) 


Initial gap g 
(nm) 


Young's 
modulus £(GPa) 


3000 


30 


100 


1000 



derived from ^ by setting the time derivative terms and 
dynamic forcing term to zero. This ordinary differential 
equation falls in the class of two point boundary value 
problems. The problem has been solved numerically in 
MATLAB environment using ROM technique and the 
solution has been verified independently by solving the 
two point boundary value problem (BVP) by collocation 
method using bvp4c solver |19j . 

For developing the ROM, the solution of © can be 
assumed as [3, Hfj ] 



u(x) 
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di<j>i(x) 
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mode shape of a straight beam under axial load N and a; 
is unknown scalar constant. The modes are orthogonal 
in nature and have been calculated analytically [20| . 

The modes were normalised such that f Q cj)fdx= 1. The 
ROM has been obtained by substituting the assumed 
solution §3§ of u{x) in ([2]). multiplying it by <p n and 
integrating from x = to x = 1. The time- independent 
form of ROM for static analysis is 



y z, j. k — 1 
l / m 

Vi=l 



1 1 ^ 

E aiajflfc / {4>i4>j)dx - J {4>l4> n )dx 



idx, for n = 1, 2...m- 



This is a system of nonlinear algebraic equation, where 
the number of equations in the system is equal to number 
of modes that has been used to develop the ROM. The 
problem has been solved with five ROM models (based 
on one to five modes). 

For construction of solution from BVP, the definite 
integral term (refer @) has been considered as a con- 
stant. The bvp4c solver requires a good initial guess of 
the solution since the solver searches the solution near 
the provided guess. Convergence problem was encoun- 
tered during solving the BVP problem for larger applied 
DC voltage near pull-in. This issue has been resolved 
using the single-mode ROM solution as an initial guess. 

Figure [2] shows the variation of normalised deflection of 
the mid-point u m id of the beam with DC voltage Vdc f° r 
five different cases of axial load N: — 0.90iVi,i, —O.BONbi, 
0, 0.50iV 6 i, 0.9(W W . Here N bl (j\r fel = A-k 2 EI/L 2 and 
Nf,i — An 2 ) is the first buckling load of a fixed- fixed 
Euler-Bernoulli beam in non dimensional form (20j . For 
each case of N, single-mode, three-modes and five-modes 
ROM results have been compared with the BVP solu- 
tion. The ROM and BVP results have also been com- 
pared with the result of Ke et al. |l8[ for the case N = 0. 
The results of current work are in good agreement with 
the published result. Ke et al. have obtained the solution 
using the finite difference method. The same nanobeam 
has also been analysed for the case N — by Ouakad 
and Younis and they observed that single-mode ROM 
provides sufficiently converged result. This is also evi- 
dent here for case N = 0, however as can be seen from 
Fig. [2] higher modes start to play an increasingly impor- 
tant role near the pull-in point when a high compressive 
stress exists in the beam. From Fig. [2j it can also be de- 
duced that as residual stress increases, pull-in voltage in- 
creases and pull-in displacement decreases. Pull-in volt- 
age for the case N = 0.90iVbi is 50% more than the case 
N = — 0.90iVbi while pull-in displacement for the case 
V = 0.90V b i is 18% less than the case N = -0.90N bl . 

The convergence of the solution has also been stud- 
ied throughout the beam deflection while Fig. [2] shows 
convergence of ROM solutions in term of u m id- Fig. [3] 
shows the static deflection curve of the beam for the case 
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FIG. 2. Static normalised mid-point deflection u m id with DC 
voltage Vdc for different values of residual stress; Nbi is first 
Euler buckling load (Inset: magnified portion of near pull-in 
point curve of the case N = 0.90^1) 



N = -0.90iV w near pull-in point at Vdc = 39.5V. This 
point is also marked with letter A in the inset in Fig. 
[2] The error in the estimation of u m id using single-mode 
ROM at this point is around 8.5%, while the error in the 
estimation of pull-in voltage for the case N = — 0.90V&1 
is around 3%. From Fig. [3] and Fig. [21 it can be deduced 
that three-modes ROM is suitable for accurate analysis 
even near pull-in. In next section, analysis of free vi- 
bration of the nanobeam at static deflected position is 
presented. 
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FIG. 3. Curvature of the deflected beam of the case N 
-0.90iV M at Vdc = 39.5V 
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IV. FREE VIBRATION ANALYSIS 

The nonlinearities in the dynamic model ([!]) cause the 
resonating frequency of the nanobeam to vary with static 
deflection under action of applied DC voltage. A linear 
vibration problem at the deflected position is required to 
be solved to determine resonating frequency. For devel- 
oping linear vibration model, the solution of (J2J can be 
assumed as 

u(x,t) = u s (x) + u d (x,t)- (5) 

The electrostatic force due to applied DC voltage de- 
termines the static displacement u s (x), whereas Ud{x,t) 
is the linear vibration component. It is assumed that 
the Ud(x,t) has very small magnitude with respect to 
static deflection u s (x) and the magnitude of Vac is a l so 
very small compared to Vdc- After substituting |S| in 
([2]) and neglecting the terms containing higher order or 
combination of Ud(x,t), derivatives of Ud(x,t) and Vac] 
a linear integro-partial-differential equation has been ob- 
tained for small vibration about the deflected position. 
The method of separation of variable has been applied by 
assuming the solution as Ud(x,t) = Udt{t)4>d{i) to obtain 
two uncoupled ordinary differential equation. The equa- 
tion in spatial coordinate is the characteristic equation, 
which provides mode shapes and natural frequencies of 
free vibration about the deflected position. The charac- 
teristic equation is 




-a 2 V£ c f' (u s )<f>d - Ud 2 (t>d = 0- 

(6) 

Here / is partial derivative of / with respect to u. This is 
an eigenvalue problem whose solution gives mode shapes 
4>d and frequencies u>a- This eigenvalue problem has been 
solved by ROM and collocation methods similar as was 
done in the static analysis. 

The solution of © can be assumed as linear combi- 
nation of modes of straight beam for developing reduced 
order model (ROM) as 

m 

4>d = 22bi<pi. (7) 

i=i 

The ROM has been obtained by substituting assumed so- 
lution ijT j) in then multiplying it by <p n and integrating 
from x = to x = 1; the resulting ROM is obtained as 

k« w n - J2 b i (Ot2V£ c J f (U s ) <f>i(f> n dx) - 
i=l \ / 

TO / 1 , 2 1 „ 1 , , 1 „ \ 

h I a.\ f u s dx ■ J <j) i 4> n dx + 2a\ J u^dx- f u s <\> n dx 

i=l V / 

= b n u>d for n = 1, 2... m- 

(8) 

This is a linear algebraic eigenvalue problem. The nat- 
ural frequencies are square roots of eigenvalues and the 



mode shapes at deflected position can be calculated from 
eigenvectors. The eigenvalue problem ^ has also been 
solved by collocation technique for verifying the results. 
Mode shapes of a straight beam have been used as initial 
guesses in bvp4c solver. 

For different values of residual stress, Fig. @]shows vari- 
ation of first resonating frequency uj\ with applied DC 
voltage Vdc- For N = Q, the results of current work has 
been compared with Ouakad and Younis result. There 
is a good agreement between results of current work and 
the published result. The observation of static analy- 
sis also holds here that the converged ROM solution can 
be obtained from three-modes ROM, and higher modes 
start to contribute in the solution near pull-in point when 
large magnitude of compressive residual stress is present 
in the beam. For the case N = —0.90^1 at Vdc = 36V, 
the error in estimation of tui from single-mode ROM is 
around 6%. 

The variation of wi versus Vdc is exploited for electro- 
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FIG. 4. Resonating frequency o?i variation with applied DC 
voltage Vdc for different values of residual stress; JV(,i is first 
Euler buckling load 

static tuning to obtain desired resonating frequency dur- 
ing operation of a NEMS resonator [6j,l2l|. Fig. 0] shows 
that presence of compressive residual stress significantly 
increases the tunability of the nanobeam resonator. In 
the case N — — 0.90iV bl , resonating frequency ui is 7.19 
at Vdc = and it increases to 24.75 at Vdc = 55V. 
Here u>i can be increased to more than three times. 

Variation of resonating frequency with residual stress 
has also been studied. Figure [5] shows this variation for 
different cases of applied DC voltage. The results are 
obtained with three-modes ROM method. The resonat- 
ing frequency u>i shows very high sensitivity on residual 
stress in absence of applied Vdc but sensitivity decreases 
under action of Vdc ■ Many resonator devices have a re- 
quirement of minimal resonating frequenc y v ariation in 
practical operating range of temperature |22l . [23j . One 
main reason for the shift of resonating frequency under 
temperature variation is build up of internal stresses due 
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to different thermal expansion coeeficients of beam and 
substrate [22[. The internal stress is similar in nature 
as residual stress. The nanobeam resonator under inves- 
tigation shows excellent frequency stability or very low 
sensitivity to internal stress in the case Vdc — 40 V. Vari- 
ation of U)\ is around 1.5% in the window of internal ax- 
ial load ±0.10-/Vbi; this is six times better than the case 
Vdc — 0. In summary, residual stress and biased volt- 
age Vdc significantly vary the resonating frequency of a 
nanobeam resonator and interplay between them may be 
utilised to increase the resonating frequency tunability 
(Fig. 2]) and to achieve frequency stability under tem- 
perature variation (Fig. [5|- The next section concludes 
the work presented in this paper. 
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based devices. In this work, a cylindrical nanobeam is 
analysed to study its static and dynamic characteristics. 
Single-mode reduced order model (ROM) is capable of 
providing accurate results when the nanobeam is free 
from compressive residual stress. However the analysis 
shows that at least three symmetric modes in ROM is 
necessary for accurate analysis near pull-in instability, 
when high magnitude of compressive residual stress is 
present. 

Residual stress has significant effects on the pull-in 
parameters and resonating frequency tunability char- 
acteristics. Presence of compressive residual stress 
decreases the pull-in voltage and increases the pull-in 
displacement. And tunability of resonating frequency 
improves significantly in presence of high magnitude of 
compressive residual stress. 

An inherent resonating frequency stability mechanism 
under temperature variation in nanobeam resonators is 
explored in this investigation. Variation in operating 
temperature causes internal stress and it is similar in 
nature to the residual stress. By studying the variation 
of resonating frequency with internal stress under action 
of constant DC voltage, it has been found that the 
resonating frequency sensitivity with internal stress 
decreases significantly at a particular DC voltage. This 
is an important property which may be considered 
during design of nanobeam resonators. 



FIG. 5. Resonating frequency uj\ variation with normalised 
axial load N/Nti for different values of applied DC voltage 
Vdc 



V. CONCLUSION 

Reduced order modelling is a simple method for analy- 
sis of electrostatically actuated nanotubes and nanowires 



ACKNOWLEDGEMENTS 

The authors gratefully acknowledge financial support 
from the Industrial Research and Consultancy Centre, 
India Institute of Technology Bombay and the Depart- 
ment of Science and Technology, Government of India. 
We are also thankful to Mandar M. Desmukh and Hari 
S. Solanki of the Tata Institute of Fundamental Research, 
India for useful discussions. 



[1] Y. T. Yang, C. Callegari, X. L. Feng, K. L. Ekinci, and 
M. L. Roukes. Zeptogram-scale nanomechanical mass 
sensing. Nano Letters, 6(4):583-586, 2006. 

[2] C. Stampfer, A. Jungen, R. Linderman, D. Obergfell, 
S. Roth, and C. Hierold. Nano-electromechanical dis- 
placement sensing based on single-walled carbon nan- 
otubes. Nano Letters, 6(7):1449-1453, 2006. 

[3] A. Subramanian, A. R. Alt, L. Dong, B. E. Kratochvil, 
C. R. Bolognesi, and B. J. Nelson. Electrostatic actu- 
ation and electromechanical switching behavior of one- 
dimensional nanostructures. ACS Nano, 3(10):2953- 
2964, 2009. 



[4] C. Rutherglen, D. Jain, and P. Burke. Nanotube electron- 
ics for radiofrequency applications. Nature Nanotechnol- 
ogy, 4(12):811-819, 2009. 

[5] M. Dequesnes, Z. Tang, and N. R. Aluru. Static and 
dynamic analysis of carbon nanotube-based switches. 
Journal of Engineering Materials and Technology- 
Transactions of the ASME, 126(3) :230-237, 2004. 

[6] H. S. Solanki, S. Sengupta, S. Dhara, V. Singh, S. Patil, 
R. Dhall, J. Parpia, A. Bhattacharya, and M. M. Desh- 
mukh. Tuning mechanical modes and influence of charge 
screening in nanowire resonators. Physical Review B, 
81(115459):l-7, 2010. 



6 



[7] H. M. Ouakad and M. I. Younis. Nonlinear dy- 
namics of electrically actuated carbon nanotube res- 
onators. Journal of Computational and Nonlinear Dy- 
namics, 5(011009):1-13, 2010. 

[8] R. C. Batra, M. Porfiri, and D. Spinello. Review of 
modeling electrostatically actuated microelectromechan- 
ical systems. Smart Materials and Structures, 16(6) :R23- 
R31, 2007. 

[9] M. M. Joglekar and D. N. Pawaskar. Closed-form empir- 
ical relations to predict the static pull-in parameters of 
electrostatically actuated microcantilevers having linear 
width variation. Microsystem Technology, 17(l):35-45, 
2011. 

[10] M. I. Younis, E. M. Abdel-Rahman, and A. Nayfeh. A 
reduced-order model for electrically actuated microbeam- 
based MEMS. Journal of Microelectromechanical Sys- 
tems, 12(5):672-680, 2003. 

[11] M. M. Joglekar and D. N. Pawaskar. Estimation of os- 
cillation period/switching time for electrostatically actu- 
ated microbeam type switches. International Journal of 
Mechanical Sciences, 53(2):116-125, 2011. 

[12] M. Rasekh, S. E. Khadem, and M. Tatari. Nonlinear 
behaviour of electrostatically actuated carbon nanotube- 
based devices. Journal of Physics D: Applied Physics, 
43(315301):1-10, 2010. 

[13] G. De Pasquale and A. Soma. Dynamic identification 
of electrostatically actuated MEMS in the frequency 
domain. Mechanical Systems and Signal Processing, 
24(6): 1621-1633, 2010. 

[14] E. M. Abdel-Rahman, M. I. Younis, and A. H. Nayfeh. 
Characterization of the mechanical behavior of an elec- 
trically actuated microbeam. Journal of Micromechanics 
and Microengineering, 12(6):759-766, 2002. 

[15] J. Kuang and C. Chen. Dynamic characteristics of 
shaped micro-actuators solved using the differential 



quadrature method. Journal of Micromechanics and Mi- 
croengineering, 14(4):647-655, 2004. 

[16] X. L. Jia, J. Yang, S. Kitipornchai, and C. W. Lim. Free 
vibration of geometrically nonlinear micro-switches un- 
der electrostatic and Casimir forces. Smart Materials 
and Structures, 19(115028):1-13, 2010. 

[17] A. Soma, G. De Pasquale, E. Brusa, and A. Ballestra. 
Effect of residual stress on the mechanical behaviour of 
microswitches at pull-in. Strain, 46(4):358-373, 2010. 

[18] C. Ke, H. D. Espinosa, and N. Pugno. Numerical anal- 
ysis of nanotube based nems devices - part (ii): Role 
of finite kinematics, stretching and charge concentra- 
tions. Journal of Applied Mechanics- Transactions of the 
ASME, 72(5):726-731, 2005. 

[19] L. F. Shampine, I. Galadwell, and S. Thomson. Solving 
ODEs with MATLAB. Cambridge University Press, 2003. 

[20] A. Bokaian. Natural frequencies of beams under com- 
pressive axial loads. Journal of Sound and Vibration, 
126(l):49-65, 1988. 

[21] B. Peng, L. Ding, and Z. Guo. Resonant modelling of 
two types of tunable carbon nanotube electromechanical 
oscillators. Micro and Nano Letters, 5(6):365-369, 2010. 

[22] M. Hopcroft, R. Melamud, R. N. Candler, W. Park, 
B. Kim, G. Yama, A. Partridge, M. Lutz, and T. W. 
Kenny. Active temperature compensation for microma- 
chined resonators. In Solid-State Sensor, Actuator and 
Microsystems Workshop, pages 364-367, Hilton Head Is- 
land, South Carolina, 2004. 

[23] R. Melamud, B. Kim, S. A. Chandorkar, M. A. Hopcroft, 
M. Agarwal, C. M. Jha, and T. W. Kenny. Temperature- 
compensated high-stability silicon resonators. Applied 
Physics Letters, 90(244107): 1-3, 2007. 



